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The spectrum of electronic states at 30000-33000 cm 1 in hafnium fluoride has been studied using 
(1+1) resonance-enhanced multi-photon ionization (REMPI) and (1+1') REMPI. Six fi' = 3/2 and 
ten iIi/2 vibronic bands have been characterized. We report the molecular constants for these 
bands and estimate the electronic energies of the excited states using a correction derived from 
the observed isotope shifts. When either of two closely spaced IR/2 electronic states is used as 
an intermediate state to access autoionizing Rydberg levels, qualitatively distinct autoionization 
spectra are observed. The intermediate state-specificity of the autoionization spectra bodes well for 
the possibility of using a selected IR/2 state as an intermediate state to create ionic HfF + in various 
selected quantum states, an important requirement for our electron electric dipole moment (eEDM) 
search in HfF + . 
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I. INTRODUCTION 



The HfF + molecular ion is identified as a promising 
candidate for an electron electric dipole moment (eEDM) 
search to test fundamental symmetries and physics be- 
yond the Standard Model The ground state 
of HfF+ is X x £+ whereas the low-lying (T e = 
977 cm -1 [Eil) metastable 3 A! state has been predicted 
to offer high sensitivity for the eEDM measurement. The 
preparation of HfF + ions in a single ro-vibrational level 
of the 3 Ai state is therefore highly desired for an eEDM 
search. 

Autoionization of Rydberg states of HfF, prepared us- 
ing the optical-optical double resonance (OODR) tech- 
nique, has been proposed as a viable way of creat- 
ing HfF + ions in a single quantum state, and has 
been demonstrated for ion-creation in the ground 
state 0- Since the ionization threshold of HfF lies at 
59462(2) cm" 1 0, 0], two ultraviolet (UV) photons of 
different frequencies are usually employed to access the 
necessary Rydberg states via the OODR technique. The 
final state of the autoionized product depends strongly on 
the intermediate state accessed by the first UV photon. 
Hence, there is a need to characterize and understand the 
neutral HfF states that lie in the vicinity of and above 
30000 cm -1 , some of which may hold promise as inter- 
mediate states for HfF + ion-creation in the desired 3 Ai 
state. 

The HfF spectroscopy published to date includes tran- 
sitions to excited states as high as 28600 cm -1 , de- 
tected by laser- induced fluorescence [8j,l9|, and resonance- 
enhanced multi-photon ionization (REMPI) [j| . This pa- 
per presents HfF spectroscopy performed in the 30000- 
33000 cm" 1 range, using both 1+1 REMPI and 1+1' 
REMPI, which complements the previously published 
spectroscopic results. 



II. EXPERIMENT 

Figure [T] shows a schematic of the experimental setup 
used to ionize a HfF molecular beam via (1+1) REMPI. 
As detailed in Ref. 0], cold (~ 10 K rotational tempera- 
ture) HfF molecules are created by ablating a Hf target 
in the presence of 1% SF 6 and 99% Ne (690 kPa back- 
ing pressure), followed by supersonic expansion through 
a 800 /im diameter nozzle. The molecular beam is col- 
limated by two 1 mm diameter skimmers separated by 
29 cm before passing through a home-built time-of-flight 
mass spectrometer (TOFMS), which is a series of six 
disks arranged along a common axis [T(|. A UV laser, 
tuned such that two photons are required to ionize HfF, 
intersects the molecular beam perpendicularly in the 
TOFMS region A. The resultant HfF + ions drift under an 
electric field of 1 V/cm toward region B. Regions B and C 
of the TOFMS are operated in the Wiley-McLaren mode 
pl| : when in region B, the ions are spatially focused by 
a small electric field of 0.25 V/cm until they reach region 
C, in which the ions are velocity- focused, by a transient 
electric field of 1 kV/cm for 500 ns, toward a microchan- 
nel plate assembly (MCP) positioned 56 cm away. The 
ion signal is enhanced by a transimpedance amplifier and 
read out by an oscilloscope. The TOFMS has a fractional 
mass resolution of at least 1/200, allowing us to resolve 
individual isotopologues of HfF + that differ by 1 amu 
(natural relative abundances in parentheses): 176 Hf 19 F 
(5.2%), m Hf 19 F (18.6%), 178 Hf 19 F (27.3%), 179 Hf 19 F 
(13.6%) and 180 Hf 19 F (35.1%). Figure^a) gives a picto- 
rial summary of all the oscilloscope traces recorded over 
multiple vibronic bands, while Fig. Oh) shows a sample 
of an oscilloscope trace for a particular wavelength of the 
UV laser. All five isotopologues of HfF + , well-separated 
in arrival time on the MCP, appear as laser frequency 
striations of different intensities. 

Most of the HfF spectra presented in this paper are 
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FIG. 1. (Color online.) Schematic of the resonance-enhanced multiphoton ionization (REMPI) experiment. 
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FIG. 2. (Color online.) (a) A pictorial summary of the os- 
cilloscope traces taken over multiple vibronic bands, (b) A 
sample of the oscilloscope trace shows the capability of the 
time-of-flight mass spectrometry (TOFMS) apparatus to sep- 
arate the spectra arising from each of the five isotopologues 
of HfF+. (c) (1+1) REMPI spectrum of 178 Hf 19 F. 



recorded using (1+1) REMPI, in which the spectral res- 
olution is limited to 0.1 cm -1 (FWHM), the linewidth 
of the UV laser. The UV photoionization radiation 
comes from a 532 nm-pumped dye laser operating with 
DCM dye, frequency doubled in a KDP crystal (100- 
200 /xJ/pulse, 10 ns, 0.1 cm^ 1 FWHM). Its wavelength 
is continuously monitored by a wavemeter calibrated 
against Ne spectral lines. 

Two of the HfF vibronic bands are recorded with high 
spectral resolution (0.003 cm" 1 ) using (1+1') REMPI, 
where two co-propagating laser pulses simultaneously 
intersect the molecular beam at the same intersection 
point. The first photon (3 /LtJ/pulse, 10 ns, 150 MHz 
FWHM) is the frequency doubled output of a home-built 
two-stage Rhodamine 101 dye cell amplifier. The dye 
cell amplifier is seeded by a continuous-wave ring dye 
laser operating with Rhodamine 610 Chloride dye. It 
is pumped by the second harmonic of a Nd:YAG laser. 
The seed laser frequency is monitored by a high-precision 
wavemeter that is regularly calibrated against the 87 Rb 
D 2 transition. The second photon in (1+1') REMPI, 
held at a fixed wavelength of 351.5 nm, is provided by 
the same dye laser as that used in the (1+1) REMPI 
setup, but operating with LDS722 dye. 



In this paper, transitions between energy levels of HfF 
are labeled as fvn/10 }, where Pq is the vibronic band 
origin in cm -1 |12j. 



III. RESULTS 

A survey scan of HfF transitions was conducted over 
the frequency range 30000-33000 cm -1 (Fig.[3Ja)), which 
is the range encompassed by the frequency doubled out- 
put of the DCM dye. Most of the strong transitions were 
scanned in detail and were found to belong to either 
an fl' = 3/2 <- fl" = 3/2 transition (Fig. E^b)) or a 



n 



1/2 



fl" = 3/2 transition (Fig. 13(c)). 



The REMPI spectra were modeled by the following 
Hamiltonian for the fl — 3/2 states (l3| : 



E{J) = T e + G l/ + BJ{J +1). 



(1) 



where T e is the electronic energy (difference in energy be- 
tween the extrapolated minima of the X 2 A 3 / 2 and elec- 
tronically excited states) and G v is the energy of the 
vibrational level v. For the Ii\/2 states, there is strong 
A-doubling and transitions from the A-doubled V 2 A 3 / 2 
ground state into both A-doublets of a given rotational 
level can be resolved, even in the low-resolution spectra. 
Since no prior information is known about the parity of a 
given member of a A-doublet, the doublets are assigned 
'aY'ri' instead of the usual spectroscopic notation 'e'/'/' 
[14| . Their energy levels can be modeled by the following 
polynomial fl5j ]: 



E a ' b (J) = T e 



7+ 



G v + BJ(J 
(-!)<' 



1) 

(P + 2?) 



(2) 



Following Ref. we assume that the 2 n spin-orbit split- 
ting is large enough for us to consider only the diagonal 

matrix clement ^Iljya |?^| 2 IIi /2^ in our description of 
the A-doublets. Among the H\/2 states, the 4 n 1 / 2 states 
can borrow oscillator strength from the 2 Hi/2 states via 
spin-orbit coupling. Since transitions into the TL1/2 up- 
per levels would look very similar to those into 2 n 1 / 2 , 
we have left out the superscript in the term symbol for 
the H1/2 levels. Based on the A-doubling, none of the 
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observed spectra belong to excited states of £ charac- 
ter, which have been previously reported in the ranges 
13800-14400 cm- 1 [9] and 19700-20000 cm" 1 Q. For 
both Hi/2 and ft' — 3/2 types of excited electronic states, 
the centrifugal distortion term, D J 2 ( J+ 1) 2 , is neglected 
because the supersonic HfF molecular beam is too cold to 
populate rotational levels beyond J = 21/2, which would 
be required for an accurate determination of D. 

Although all five isotopologues were observed, only the 
two most abundant ( 178 Hf 19 F and 180 Hf 19 F) were an- 
alyzed. (In Section IIVD| we extend for selected bands 
our analysis to the 177 Hf 19 F and 179 Hf 19 F isotopologues.) 
For a given transition, the two isotopologues' REMPI 
spectra are simultaneously fit to a contour described by 
a common set of fit parameters: temperature, intensity, 
band origin for 180 Hf 19 F ( 18 % 



180 



t: 



_180ji// _j_180 qi 



180 G'l„), band origin for 178 Hf 19 F ( 178 ^ ), lower state 



rotational constant {B"), excited state rotational con- 
stant (£?') and, where applicable, the A-doubling con- 
stant (p + 2q)/2. Since the rotational constants of a 
given state for the two isotopologues are related to each 
other by the inverse ratio of their reduced masses \i (e.g., 
18Q B/ 17S B = 178 m/ 18 V), the B-value ratios are held 
fixed for the simultaneous contour fit. The resolution 
of our data is insufficient to detect a difference between 
the (p + 2q) /2 terms for different isotopologues. 

The molecular constants for the observed HfF tran- 
sitions are summarized in Table Mil The numbers in 
parentheses indicate the ltr-standard deviation in the last 
digit. The standard deviation for each fit parameter was 
calculated from a histogram of fit parameters obtained by 
bootstrapping the residuals of the fit [l6[ . The bootstrap- 
ping process adds the residuals, obtained after a single 
fit iteration, randomly to the fitted data to form a new 
data set, which is then re-fit against the model; this pro- 
cess is repeated a hundred times to get a hundred sets 
of fit parameters for the histogram. For the bands ex- 
amined in high-resolution (1+1') REMPI scans, (1+1) 
REMPI data for the same bands was also available for 
comparison. We find good agreement between the values 
of vq from the coarse (1+1) REMPI spectra and from 
the (1+1') REMPI spectra, where the uncertainty in v$ 
is the specified accuracy of the wavemeter, i.e. 0.1 cm -1 . 
The global systematic error of 0.1 cm -1 is determined 
by how often the wavemeter is calibrated against the Ne 
spectral lines. For each vibronic band, because the spec- 
tra of all isotopologues are simultaneously recorded, the 
isotope shifts from the (1+1) REMPI scans and (1+1') 
REMPI scans agree well when the error calculated by the 
bootstrap method is at least 0.01 cm -1 . Similarly, the 
uncertainties in B' and B" appear to be dominated by 
the bootstrap statistics when they are at least as large as 
0.001 cm -1 . When the purely statistical standard devia- 
tions in the isotope shift or {£?', B"} fall below 0.01 cm -1 
and 0.001 cm -1 respectively, the uncertainties reported 
in Table IIIII have been increased to reflect an estimate 
of the corresponding systematic errors. We find that the 
purely statistical uncertainty estimated using the boot- 



strap method accounts well for the uncertainty in the 
A-doubling parameter. 



IV. DISCUSSION 
A. Lower states of observed transitions 

From previous work, the electronic ground state of 
HfF is 2 A 3 / 2 and for v — 0, the rotational constants 
are B =0.284 001(7) cm" 1 and 0.283 668(6) cm" 1 for 
the 178 Hf 19 F and 180 Hf 19 F isotopologues respectively H. 
From our own laser-induced fluorescence studies [9j, we 
get the isotope-averaged rotational constants to be £?„ = o 
= 0.2836(4) cm -1 , B v=1 = 0.2822(10) cm" 1 and B v=2 
= 0.2791(12) cm" 1 . We found 9] that in our appara- 
tus a majority of the HfF molecules are created in the 
ground vibronic state, although populations in multiple 
lower vibrational levels can be retained, athermally, in 
the supersonic expansion. 

In this paper, all characterized transitions appear to 
originate from the 2 A 3 / 2 ground state of HfF. In prin- 
ciple, the vibrational assignment of the lower state can 
be inferred from a precise fit to the lower state's rota- 
tional constant and comparison against the values ob- 
tained in Ref. [9]. The vibrational quantum number 
of the lower state was indeed obtained from our high- 
resolution (1+1') REMPI data (Fig.^d)), which shows 
that some bands we observe come from the ground vi- 
bronic state. However, the (1+1) REMPI data are too 
noisy, too broad in linewidth, and have population in too 
few rotational states to permit a precise enough determi- 
nation of B". 



B. A tool for understanding vibronic bands 

The spectra of HfF in the vicinity of 30000 cm -1 are 
complicated; tools are needed to sort out the spectra and 
to group together excited states that belong to the same 
electronic orgin. For most of the vibronic bands, we did 
not have available to us sufficiently precise determina- 
tions of both B' and B" to determine the vibrational 
numbering of upper and lower state levels. Instead, we 
used the isotope shift, Au = 180 ^o — ™ v 0i to determine 
an estimate of the electronic energy, T' e , for each excited 
vibrational level. 

Assuming that the electronic and vibrational energies 
are separable, we can write the vibronic band origin for 
the isotopologue *Hf 19 F as 



(iW _ (i) T n , {i)r<i _ (i)n" 



(3) 



and, assuming an approximately harmonic potential, we 
can write 
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FIG. 3. (Color online.) (a) Coarse (1+1) REMPI scan of excited HfF states near 30000 cm -1 . The ion signal shown is the 
integrated ion signal for all five HfF isotopologues. Two main types of transitions have been identified in this survey scan: 
Si' — 3/2 4— SI" = 3/2 (marked by asterisks, typical spectrum shown in (b)) and II1/2 4- SI" = 3/2 (marked by arrows, typical 
spectrum shown in (c)). (b) Detailed (1+1) REMPI scan of the {31.89} vibronic band. The red smooth curve is a contour fit 
to the data shown as connected black dots, (c) Detailed (1+1) REMPI scan of the {32.98} vibronic band. The blue smooth 
curve is a contour fit to the data shown as connected black dots, (d) High resolution (1+1') REMPI spectrum of the {32.32} 
vibronic band. Individual upper state rotational levels of opposite parity can be resolved. 



TABLE I. Summary of molecular constants (all in cm ) for observed transitions. The numbers in parentheses denote the la 
uncertainties, assigned as explained in the text. 
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0.282(3) 


0.249(3) 
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-1.838(10) 
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where M/x is the reduced mass for isotope (i). A similar 
relation applies for G"„, thus if we define 180 T" = 0, we 
combine Eqs. (|3J) and (d} to give an estimated value of 
the electronic energy, 180 T e : 

= 18 % + 77 [Av 

- ( 180 Z ~ 178 Z) + ( 180 t; - 178 T e ")] , (5) 



where 77 = y /5 ^/(v /iM M- V 17 ^ 1 ) ~ 1861 - Then > as- 
suming that the electronic contribution to the isotope 
shift is much smaller than the vibrational contribution 




180 T e ' = 18 % + ^. (6) 

Even in the case that the electronic contribution to 
the isotope shift ( 180 T^ - 178 T e ' - 180 T e " + 178 T e ") is not 
much smaller than the vibrational contribution { 1W G' V , — 
178 G' U , - 180 G"„ + 17S G'l„), two bands that share com- 
mon upper and lower electronic states should get, from 
Eq. [JJJ very similar values for 180 Tg, even if v" and v' are 
different for the two bands. Again, this depends on the 
assumption of harmonic potentials and of separable elec- 
tronic and vibrational contributions to the band energy. 
Figure U] is a graphical summary of the observed SI' = 
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3/2 <- X 2 A 3/2 bands and n 1/2 <- X 2 A 3/2 bands in HfF. 
Part (a) of the figure shows the isotope shifts, where an 
isotope shift is estimated to be about -0.3 cm -1 for a 
change in vibrational quanta of v 1 — v" = +1. Some 
isotope shifts of anomalously large magnitude (1.5 cm -1 
to 3 cm -1 ) are further discussed in Section TlVDI Part 
(b) shows the estimated electronic energy versus band 
origin, in which transitions that share the same pair of 
lower and upper electronic states are expected to line up 
horizontally. In Fig. 0Jb) we see two bands {31.89} and 
{32.47} that, while their band origins are separated by 
579 cm -1 , have 180 Tg values that differ by only 13 cm -1 
(identical within measurement uncertainty). It is likely 
that these two bands share a common electronic transi- 
tion with 180 T e ' w 31840 cm -1 , as indicated by the y-axis. 
We note that 579 cm -1 would be a reasonable excited 
state vibrational spacing. We can see, from the spacing 
of data points in the y-direction in Fig. QJb), that there 
are at least four distinct f2' = 3/2 electronic levels: near 
27200, 28150, 31850 and 32400 cm" 1 . As a caveat, we 
note that several of the bands have unreasonably large 
isotope shifts. We discuss possible causes for anomalous 
isotope shifts in Section TlV Dl below. 



For the H\/ 2 ^— X 2 A 3 / 2 transitions, from Fig. SJb) we 
can see immediately that there are multiple distinct Tli/ 2 

electronic levels: one with 180 T e ' » 27500 cm -1 , at least 
one with 180 Tg « 29600 cm -1 , and multiple levels with 
180 Tg in the range 31500 ± 300 cm -1 . One cannot say for 
certain that the three bands with 180 Tg w 29600 cm" 1 
share a common electronic transition. With reference 
to Table I, we see the {32.34} band has a A-doubling 
constant quite distinct from that of the {32.91} band, 
which suggests either two distinct electronic transitions 
or, more likely, a local perturbation of one of the bands. 
For the cluster of six O' = 1/2 bands with 180 T^ be- 
tween 31200 and 31800 cm -1 , the spread in 180 5^ is too 
large for the bands to share a common electronic tran- 
sition. Moreover, the respective rotational constants for 
the bands {32.22} and {32.32} are determined so pre- 
cisely that we can say with some confidence that these 
two bands share a common lower level vibrational level 
[v" = 0), but belong to distinct upper vibrational levels. 
Since the values of vq are too close together to make it 
likely that the upper levels are adjacent vibrational levels 
of the same electronic state, it seems quite likely that the 
two bands terminate in distinct electronic upper levels. 
For the bands {32.19} and {32.86}, the vibronic band ori- 
gins differ by 668.41(14) cm -1 , which is within 2a of the 
1-0 vibrational energy splitting of the X 2 A 3 / 2 ground 
state @. Further, both bands fit to the same B' whereas 
one band fits to a value of B" that is consistently smaller 
than that of the other band. This strongly suggests that 
the bands {32.19} and {32.86} arise from the v" = 1 and 
v" = ground vibrational levels, respectively. 



C. Rotational line strengths 

In the (1+1) REMPI experiments, the transition ex- 
cited by the second photon is typically left unsaturated, 
so the observed REMPI line strengths are actually a con- 
volution of the line strengths for both the first and second 
step in REMPI. We have observed in other work [7j that 
the above-threshold ionization spectrum is highly struc- 
tured. This structure corresponds to z-'Ryd > levels 
that autoionize into v + = vibrational levels of HfF + . 
Qualitatively, the Honl-London expressions for Hund's 
case (a) molecules work well to describe the observed ro- 
tational line intensities. This is the case for all observed 
bands except two tt' = 3/2 <- X 2 A 3/2 bands: {30.73} 
and {32.47}. 

To fit the line intensities in a given vibronic band, we 
included terms that described possible interference ef- 
fects between parallel — ft" = 0) and perpendicular 
(O' -O" = ±1) transitions For example, a nominal 
O' = 3/2 (for clarity, also denoted here as £l' N ) excited 
state could also possess some Q' = 1/2 (= fl' N — l) and/or 
D,' = 5/2 (= Q' N + 1) character from mixing (via the 
Hamiltonian term -SJ^L*) with nearby "dark" states. 
When the excited state is no longer a good Hund's case 
(a) state due to these admixtures, transitions between 
the 2 A 3 j 2 ground state and all three admixed fi' compo- 
nents are allowed. The line strength S(J',Cl' N ; J"Cl") for 
a given rotational line is given by: 



s(j',n' N ; j"n") oc 
j' i 



Mo 

+/i 7 



J" 



n' N (n' N - ci") ci" 



j' i 
-(o^ + i) (Q^ + i-n") 

j' i 
-(o^-i) (n^-i-n") 



j" 

Q" 

j" 
n" 



(7) 



where fi a ,fj,ff, and /x 7 refer to the transition dipole ma- 
trix elements between the ground state and the admixed 
excited states (of Q' N , fl' N + 1, and fl' N — 1 character, 
correspondingly). Note that, since the fi a , [ip, /i 7 ampli- 
tudes are summed and then squared, interference effects 
are present. These matrix elements are allowed to vary 
in the contour fits. If the second and third terms are ig- 
nored, the line stren gth expression reduces to the normal 
Honl-London factor [17] . 

For the H\/ 2 excited states, the Ail = ±1, transition 
selection rule requires the last term to be zero. Between 
the two remaining transition dipole matrix elements, \x a 
tends to dominate for all of the observed IIj/2 X 2 A 3 / 2 
bands. Similarly, [i a tends to be the dominant dipole 
matrix element for all of the f2' = 3/2 ■< — X 2 A 3 / 2 bands 
except for two, {30.28} and {32.47}. Figure [5] compares 
the contour fits performed using the Honl-London expres- 
sions versus Eq. (J7|) to describe rotational line intensities 
for each of the anomalous bands. In both cases, the in- 
corporation of interference effects in Eq. {7} gives rise to 
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a better fit. In particular, the contour fits indicate an 
admixture from a nearby tt' = 1/2 state (where /i 7 is 
dominant), which suggests that the observed upper state 
is of nominal 2 n 3 / 2 character, the transition "brightness" 
from the ground state has been lent to it from a nearby 
H1/2 state. In fact, the rotational line strengths imply 
that at least one of the two bands ({32.47}) is located 
very close to a II1/2 <~ X 2 A 3 / 2 band, as seen in Fig. @] 
The observed admixture of electronic states is the reason 
for using Hund's case (c) notation to specify the £1' = 3/2 
excited states. 



D. Large isotope shifts 

As shown in part (a) of Fig. 2J several bands were 
observed to have isotope shifts more negative than - 
1.6 cm , which would correspond to the upper state of 
the transition having vibrational energy that exceeds the 
lower state vibrational energy by more than 3000 cm -1 . 
If we assume a typical upper-state vibrational constant 
of ~ 600 cm -1 , a 3000 cm -1 vibrational change im- 
plied a change in vibrational quanta, Av ~ 5. Deriving 
bond lengths from observed rotational constants, one cal- 
culates extremely small Franck-Condon factors for such 
transitions. This poses a mystery: how were we able to 
excite these transitions? One possible explanation could 
be that interactions between potential curves, say an an- 
ticrossing of some sort, generate a few levels described by 
anomalously large vibrational constants, as depicted in 
Fig. [H[a). An upper-state vibrational constant as large 
as 1000 cm -1 , for instance, could yield an acceptably 
large Franck-Condon overlap for the corresponding Av = 
3 transition. But it seems unlikely there would be multi- 
ple electronic states with such exotically large curvatures 
in their potential curves. Moreover, in Fig. [4j the x-axis 
spacing between some pairs of points that are clumped 
along the y-axis suggests that at least some of the large- 
isotope-shift bands have relatively modest upper-state vi- 
brational spacing, perhaps 450-600 cm -1 . 

As an alternative hypothesis to the picture of tran- 
sitions involving large changes in the number of vibra- 
tional quanta, the anomalous isotope shifts could arise 
from isotope-specific accidental degeneracies between two 
mutually perturbing excited states, each of which have, 
in the absence of perturbation, smaller isotope shifts 
(Fig. EKb)). In this case, the perturbing state would not 
only cause the 178 Hf 19 F and 180 Hf 19 F states to split apart 
from each other, but also cause all of the other isotopo- 
logues to follow an irregular splitting that is not linear 
with ^/l/W/x. Figure [7] displays the isotope splitting be- 
tween the four most abundant isotopologues for several 
bands with large isotope shifts. The isotopologues are 
found to follow a linear energy spacing relative to the in- 
verse square root of their reduced masses, which is what 
we expect in the absence of such an isotopologue-specific 
perturbation and where the observed isotope shift is de- 
termined primarily by ylj W [i vibrational energy shifts. 



The high visibility of bands that had large isotope 
shifts could also be due to spin-orbit interactions be- 
tween vibrational levels that belong to different excited 
states, as shown in Fig.[5tc). The observed excited state 
of a high vibrational quantum number, v' d , could have 
been an initially "dark" electronic state that became ob- 
servable by spin-orbit interaction with a "bright" elec- 
tronic state of much lower vibrational quantum number, 
v' h . This would provide the latter with decent Franck- 
Condon overlap with the lower state. Such a spin-orbit 
interaction would be proportional to the vibrational over- 
lap matrix element, {v' d \v'^) [H|], which is estimated to be 
large only for v' d w i/' b ± 1 for all the observed excited 
states in HfF. 

One plausible explanation for the large isotope shifts 
is that there may be a third state that induces an indi- 
rect interaction between the observed state and another 
"bright" state, where the "bright" state is of a signifi- 
cantly lower vibrational quantum number and has large 
Franck-Condon overlap with the lower state (Fig. HJd)). 
It is likely that there are other plausible explanations for 
the observation of such anomalously large isotope shifts, 
and we invite the interested spectroscopist to offer his or 
her ideas. 



E. Intermediate states for OODR autoionization 

Among the six U 1/2 <- X 2 A 3/2 bands with 180 T^ in 
the range 31200-31900 cm" 1 , two of the bands ({32.22} 
and {32.32}) have been scanned at high resolution, using 
the (1+1') REMPI technique. Their distinct A-doubling 
constants, coupled with the close proximity of the elec- 
tronic energy levels, indicate that they must belong to 
electronically distinct states in order to exhibit such small 
repulsion. Their dominant electronic configurations must 
differ by at least two spin-orbitals, e.g. sdS( 3 Ai)n£X ver- 
sus s 2 ( 1 E+)n'£ / A', where sdS^A^ or s 2 ( 1 E+) refers to 
the core configurations, and n£X or n'l'X' refers to an ad- 
ditional electron in a more highly excited orbital. These 
two configurationally distinct states could potentially be 
used as intermediate states to access different states of 
HfF + when performing OODR autoionization, which is 
relevant to our goal of preferential population of the 
metastable 3 Ai state rather than X X Y, + for the eEDM 
experiment. These two vibronic states have indeed been 
observed to yield very different autoionization spectra, as 
shown in Fig. [5J We are in the process of characterizing 
the electronic states of the ions so produced. 

V. CONCLUSIONS 

A plethora of HfF bands in the 30000-33000 cm" 1 re- 
gion have been observed using (1+1)REMPI and (1+1') 
REMPI. We have characterized six Q' = 3/2 «- X 2 A 3/2 
and ten II1/2 <— X 2 A 3 / 2 vibronic bands. To sort out 
the spectra, we used the isotope shift for a given band 
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to determine the electronic energy, 180 T e ', for the upper 
electronic state. This method of grouping bands only 
works for bands where the potential energy curves are 
fairly well-approximated by a harmonic potential. Two 
bands exhibit interference effects between parallel and 
perpendicular transition moments through their rota- 
tional line strengths. Several bands with anomalously 
large isotope shifts had intensities much larger than pre- 
dicted based on the expected small Franck-Condon fac- 
tors for transitions from the low-v" lower state. Among 
the six ^ — A" 2 A3/2 bands with electronic energy off- 
sets crowded in the vicinity of 31300-31800 cm -1 , there 
are at least two electronically distinct states. When such 
configurationally distinct states are used as intermediate 
states in the OODR preparation of Rydberg states, these 



states provide access to at least two possible routes for 
creating distinct HfF + electronic states after autoioniza- 
tion decay, which will be important for the selective for- 
mation of ionic HfF + in the desired 3 Ai quantum state 
for the eEDM experiment. 
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FIG. 4. (a) Isotope shifts (triangles denoting Q' = 3/2 <— 
X 2 A 3 / 2 transitions and circles denoting ITj/2 ^~ X 2 A 3 / 2 
transitions) and (b) calculated electronic energies for the 
fl' — 3/2 X 2 A 3 / 2 transitions (dashed lines terminating on 
the left axis) and H1/2 X 2 A 3 / 2 transitions (dashed lines 
terminating on the right axis). The spread in 1S0 T e ' values 
indicates at least four distinct fl' =3/2 electronic states and 
at least six distinct H1/2 electronic states. 
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FIG. 5. (Color online.) Contour fits (red) to the Q,' = 3/2 <s— X 2 A 3 / 2 {30.73} and {32.47} bands (data shown as connected 
black dots). In (a, c), the contour fits are constrained to use Honl-London expressions for the rotational line strengths, which 
are shown here to be inadequate for a good description of the observed line intensities, especially for those in the Q branch. 
In (b, d), the contour fits are allowed to take into account possible interference effects between parallel and perpendicular 
transition moments. In both cases, /i 7 dominates over the other transition dipole matrix elements, indicating that the nominal 
Q,' — 3/2 upper states in both cases appear to possess an admixture from nearby SI' = 1/2 states. 
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FIG. 6. Four hypotheses for the observation of large isotope 
shifts, (a) Avoided crossing of potential energy curves, giving 
rise to anomalously large vibrational spacings. (b) Isotope- 
specific accidental degeneracies between two mutually per- 
turbing excited states, ip[ and ip'2. In this picture, the isotope 
splittings would not be linear in -^/l/W/i. (c) Local perturba- 
tion between a high-i/^ level of a "dark" electronic state and a 
low-i^ level of a "bright" electronic state. The higher u' d level 
becomes observable by borrowing brightness from the \ow-v' b 
level, but it has the normal isotope shift of a high-z^ level, (d) 
Indirect coupling between the high-j/^ level and \aw-v' b level. 
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FIG. 7. For the vibronic bands {30.73}, {32.34}, {32.91}, and {32.98}, the isotope shifts are found to vary linearly with the 
inverse square root of the isotopologue reduced masses, which suggests that anomalously large isotope shifts are indeed due to 
a large Av and not to isotope-specific perturbations. 
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FIG. 8. OODR autoionization spectra measured by tuning the first photon to access a particular excited rotational state 
(J' = l/2(a)) using the vibronic bands (a) {32.22} and (b) {32.32} (ion signal inverted for clarity), then scanning the second 
photon to map out the spectrum of autoionizing Rydberg states accessible from that intermediate state. Each case gives a 
unique set of autoionization resonances, suggesting that the ionic core of the autoionizing Rydberg state differs between the 
two cases. 



